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ABSTRACT 

We examine the evolution of variously-sized radio galaxies [i.e., compact symmetric objects (CSOs), 
medium-size symmetric objects (MSOs), Fanaroff- Riley type II radio galaxies (FRIIs)], by comparing 
the relation between the hot spot size and the projected linear size with a coevolution model of hot 
spots and a cocoon. We take account of the deceleration effect by the cocoon head growth. We find 
that the advance speed of hot spots and lobes inevitably show the deceleration phase (CSO-MSO phase) 
and the acceleration phase (MSO-FRII phase). This is ascribed to the change of the power-law index 
of ambient density profile in the MSO phase (' 1 kpc). It is also found that the cocoon shape becomes 
nearly spherical or disrupted for MSOs, while an elongated morphology is predicted for CSOs and FRI 
Is. This seems to be consistent with the higher fraction of distorted morphology of MSOs than that of 
CSOs and FRIIs. Finally, we predict that only CSOs whose initial advance speed is higher than about 
0.1c can evolve into FRIIs, comparing the hot spot speed with the sound speed of the ambient medium. 

Subject headings: galaxies: active — galaxies: evolution — galaxies:jets — galaxies:ISM 



1. INTRODUCTION 

According to the standard model for large-scale radio 
sources (Blandford & Rees 1974), a pair of relativistic jets 
transport away the bulk kinetic energy from the central 
compact region close to a central supermassive black hole 
(SMBH) to - 100 kpc - 10^ kpc scale radio lobes. This 
bulk kinetic energy is dissipated by the strong terminal 
shocks which are identified as hot spots. Then, shocked 
relativistic plasma expands sideways and envelopes the 
whole jet system and this is a so-called cocoon (Scheuer 
1974). Now, how young AGN jets (or radio sources) evolve 
into the extended radio sources (e.g., FRII radio galaxies) 
is one of the primal problems in astrophysics. In order 
to clarify this issue, it is important to understand the na- 
ture of small, young progenitors of FR lis. The study 
of young progenitor of FRIIs is also a good opportunity to 
get information about the central regions of galaxies, espe- 
cially their mass density distribution. The coevolution of 
SMBHs and galaxies has been discussed by many authors 
who say that AGN feedback may be a key mechanism in 
regulating the growth of SMBHs and also in preventing 
further star formation of galaxies; These authors use an 
analytical approach (e.g., Silk & Rees 1998; King 2003; 
Wyithe & Loeb 2003; Murray et al. 2005) and a semi- 
analytical approach (e.g., Granato et al. 2004; Kawata 
& Gibson 2005; Di Matteo et al. 2005; Okamoto et al. 
2007). However, it has been unclear how we should treat 
AGN feedback associated with AGN jets, because we do 
not understand how small, young radio sources evolve into 
FREs. Therefore, for understanding AGN feedback, it is 
essential to explore conditions under which small, young 
radio sources can pass through a dense ambient medium 



within host galaxies. 

Recently, a large number of compact, bright double- 
lobe radio sources, so-called "compact symmetric objects 
(CSOs)" whose linear size is < 1 kpc have been discov- 
ered (Wilkinson et al. 1994; Fanti et al. 1995; Readhead 
et al. 1996a, b). Concerning the origin of CSOs, two 
scenarios were proposed. One is a so-called "frustrated 
jet scenario" in which the ambient medium is so dense 
that the jet cannot break through, so that sources are old 
and confined (van Breugel, Miley & Heckman 1984). The 
other is a "youth radio source scenario" in which CSOs are 
the young progenitor of FRII radio galaxies (e.g.. Philips 
& Mutel 1982). The former scenario looks less likely be- 
cause recent multi- wavelength studies show that the ambi- 
ent medium of CSOs is not significantly different from FRI 
I radio galaxies (O'Dea et al. 1998 and references therein). 
On the other hand, the most compelling evidence for the 
youth jet scenario is found in the recent observations of 
advance speed of hot spots of several CSOs (Owaiank, 
Conway & Polatidis 1998; Taylor et al. 2000; Tschager 
et al. 2000; Giroletti et al. 2003; Polatidis & Conway 
2003; Gugliucci et al. 2005; Nagai et al. 2006; Gugliucci 
et al. 2007; Luo et al. 2007). They showed that the sepa- 
ration velocity is typically ^ 0.1c, indicating a dynamical 
age of ~ 10^"'* yr which implies that CSOs are possible 
candidates as the progenitors of FRIIs. 

Based on the youth scenario, the evolution of radio 
sources has been considered by a number of authors (Car- 
valho 1985; Begelman & Cioffi 1989, hereafter BC89; Falle 
1991; Fanti et al. 1995; Readhead et al. 1996a; Begelman 
1996; De Young 1997; Kaiser & Alexander 1997; Alexan- 
der 2000; Snellen et al. 2000; Perucho & Martf 2002; 
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Kawakatu & Kino 2006, hereafter KK06). For simplicity, 
one has supposed a constant advance speed of hot spots, 
or a constant aspect ratio of the cocoon (i.e., a self-similar 
evolution) for a long-term evolution of radio-loud AGNs. 
Do these treatments reflect the actual evolution of radio 
sources ? The answer is not trivial. According to previ- 
ous observations, a stronger interaction with the ambient 
medium has been reported for more compact radio sources 
(e.g., Gelderman & Whittle 1994; de Vries et al. 1997; 
Axon et al. 2000; Gupta et al. 2005; Hoft et al. 2008). 
Thus, the advance speed of hot spots might decelerate via 
the strong interaction in host galaxies. On the evolution of 
cocoon morphology, recent observations have shown that 
the fraction of distorted radio morphology of medium-size 
symmetric objects (MSOs) is higher than CSOs and FRI 
Is (e.g., Saikia et al. 1995; O'Dea 1998; Dallacasa et al. 
2002a, b). Observationally, it is unclear whether the as- 
sumptions (constant hot spot velocity or constant aspect 
ratio) are reasonable. Also, the validity of these assump- 
tions is still under debate theoretically (De Young 1997; 
Komissarov & Falle 1998; O'Dea 1998; Carvalho & O'Dea 
2002a, b; Perucho & Martf 2003: hereafter PM03). There- 
fore, in order to answer the above question, it is essential 
to build up an appropriate model of radio sources with- 
out assuming the constant aspect ratio and the constant 
advance speed of hot spots presented by KK06. 

Recently, it has been found that the power law index 
for the evolution of hot spot size changes at the transi- 
tion between the interstellar medium and the intergalactic 
medium, i.e., ^ 1— lOkpc ( Jeyakumar & Saikia 2000: here- 
after JOO; PM03). Since the hot spot is one of the most 
important ingredients in the whole jet system, the evolu- 
tion of hot spot size would reflect the dynamical growth 
of radio sources of various scales. However, it was difficult 
to derive the dynamical evolution of radio sources from 
previous work (JOO and PM03) because of the lack of spa- 
tial resolution, the observational bias and small sample 
of radio sources. Therefore, we should not only compile 
a larger sample of CSOs, MSOs and FRUs, but also be 
careful about the data quality and the observational bias. 
Then, by comparing the direct comparison with KK06, we 
examine the deceleration and acceleration of hot spots. 

The paper is organized as follows: In §2, we first elu- 
cidate whether the power-law index change at '-^ 1 kpc is 
a solid result for the evolution of hot spot size (§2.1). To 
this end, we compile a larger sample of CSOs, MSOs and 
FRn sources than previous work (JOO; PM03) by consid- 
ering the observational bias and being careful about the 
data quality. We also discuss the advance speed of the 
hot spot (§2.2). In §3, to interpret these observational 
data, we briefly review the physical model of hot spot evo- 
lution (KK06) with the aid of cocoon dynamics (Kino & 
Kawakatu 2005; hereafter KK05). In §4, from the direct 
comparison with KK06 we reveal the deceleration and ac- 
celeration of advance speed of hot spots, and the evolu- 
tion of cocoon morphology. In §5, we predict the fate 
of CSOs by comparing the evolution of hot spot velocity 
with the sound speed of the ambient medium. Finally, we 
discuss the FRI/II dichotomy and AGN feedback due to 
the cocoon (i.e., AGN bubble) expansion. Section 6 is de- 
voted to our summary. Throughout this paper, we adopt 
Ho = 72 km s'^Mpc^^ and qo = 0.5. 



2. OBSERVATIONAL PROPERTIES OF HOT SPOTS 

Following the recent work of Nagai (2007), we com- 
piled variously-sized radio galaxies and examined phys- 
ical quantities of the hot spot (Table 1). FRII sources 
('ii ^ 10 kpc, where is linear size of radio source) in our 
sample are selected from FRII radio galaxies in the 3CR 
catalogue (Hardcastle et al. 1998; Fernini et al. 1993; Bri- 
dle et al. 1994; Gilbert et al. 2004; Kharb et al. 2008). 
MSOs (1 < ^ 10 kpc) are selected from the samples 
in Dallacasa et al. (2002a) and Dallacasa et al. (2002b). 
We also included Comact-Steep Spectrum (CSS) sources 
having double-lobe-like structure in Fanti et al. (1985) 
and Sanghera et al. (1995) in our MSOs sample. CSOs 
(< 1 kpc) are selected from the samples in Readhead et 
al. (1996b), Taylor et al. (2000), Wang et al. (2003), Po- 
latidis et al. (2003), Maness et al. (2004) and Nagai et al. 
(2006). In addition to these data, we measured hot spot 
parameters of a CSO (B 1943-1-546) by analyzing the data 
obtained from the NASA/IPAC Extragalactic Database 
(NED). Sources with a complex and highly distorted struc- 
ture have been excluded from the above samples. To avoid 
the difference in estimation of physical quantities between 
the hot spot and the counter hot spot due to the projection 
effect, we focus on sources with relatively symmetric radio 
lobes in this work. Finally, we compiled a total of 117 ra- 
dio sources. For the sources with multiple hot spots, we 
followed the method of hot-spot identification described in 
the original papers. If there was no information of identifi- 
cation method available in the original papers, we referred 
to the most prominent component as the hot spot. Based 
on these data, we investigate how the size of the hot spot 
changes with the distance from the core in §2.1. We also 
discuss the advance speed of the hot spot in §2.2. 




Distance from the core [kpc] 

Fig. 1. — The relation of hot-spot size (rns) and hot-spot 
distance from the core (Ih)- Crosses with arrows indicate upper 
hmit. Sohd line corresponds to the best-fit for the sources where 
10""^ kpc< /h ^ t kpc whereas broken line corresponds to that for 
the sources where 1 kpc< < 10^ kpc. Note that upper-limit data 
were not included in the fittings. 
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2.1. Size 

Figure 1 shows the size of a hot spot, rjjs with re- 
spect to Zh, together with the best Hnear fit on the log- 
log plane (logy = alogx -I- h). We estimated rns as 
(^maj + ^min)/2 where 0,naj and ^min are full width at half 
maximum (FWHM) of major and minor axes of the com- 
ponent obtained by Gaussian model fittings, respectively. 
The slope is changed around ~ 1 kpc, and then the best 
fit value of a for Zh < 1 kpc is a = 1.34 ± 0.24, while 
a = 0.44 ± 0.08 for Zh > 1 kpc. This tendency is almost 
consistent with PM03, but a large number of sources in 
our sample allow us to confirm the tendency more clearly. 
One could think the apparent slope change is caused by 
observational bias. Here we consider the possible uncer- 
tainties of the estimation of the hot-spot size and the slope 
change. 

Firstly, the sizes of hot spots of larger sources could be 
underestimated due to a lack of sensitivity for diffuse emis- 
sion. In this context, there are two major candidates for 
underestimation of size. One is the missing fiux for the 
diffuse emission due to the absence of short baseline of 
the interferometer. Generally, this effect is more serious in 
VLBI (e.g., VLB A, EVN) observation than in interferom- 
eter (e.g., VLA) observation. The hot-spot sizes of most 
sources with Zh > 1 kpc are estimated using data taken 
with VLA while that with Zh < 1 kpc are estimated using 
data taken with VLBL For the sources with Zh > 1 kpc, 
the missing flux is not significant because almost the to- 
tal extent of the radio lobe has been, in fact, detected in 
most sources in our sample. The other candidate is the 
image dynamic range. Poor dynamic range could contam- 
inate the diffuse emission so that the hot-spot sizes could 
be underestimated. The image dynamic range strongly 
depends on the distribution of sampled visibilities on the 
uv-coverage if we assume each visibility has no significant 
errors (Perley 1999). Denser and more uniform sampling 
of visibilities on the uv-coverage increases the image dy- 
namic range. Thanks to dense and uniform sampling of 
the visibility, the image dynamic range of VLA observa- 
tion attains up to several thousands (e.g., Hardcastle et 
al. 1997), which is much higher than typical VLBI obser- 
vations. Therefore, we can safely exclude the possibilities 
of underestimation of the hot-spot sizes for larger sources. 

Secondly, the angular resolution relative to the largest 
angular size of radio source (rtb) possibly affects the esti- 
mation of the hot-spot size. With high n\, value one could 
expect to see smaller-scale structures within the hot spots, 
which leads to a decrease in the sizes of hot spots. FRK 
sources in our sample have observed with nb that is a few 
tens of times higher than n\, of CSOs and MSOs in our 
sample. In this case, smaller scale structure within hot 
spots could be seen for FRU sources than those of CSOs 
and MSOs. However, if CSOs and MSOs are observed with 
the same order of rib for FRH sources, the hot-spot sizes 
of CSOs and MSOs should decrease. This effect is more 
serious in smaller sources. As a result, the slope in CSOs- 
MSOs phase would be more inverted (a > 1.34 ±0.02) but 
the starting point of slope change would be unchanged. In- 
tensive study by JSOO found a similar change in slope, but 
the change occurred around ^ 20 kpc. They selected the 
sample observed with similar nb among CSOs, MSOs, and 
FRH sources. However, their sample contains many FRII 



sources where their hot-spot sizes have been determined as 
only upper limits due to lack of angular resolution. This 
would make the starting point of the slope change unclear. 

Lastly, we comment on the projection effect. To be 
precise, our estimated size of the hot spot is not intrin- 
sic size but the size projected onto the celestial plane. 
However, our sample was restricted to the sources with 
relatively symmetric lobes. We excluded the one-sided 
sources. Thus the difference in the jet axis is not so large 
among the sources, and affects the estimation of the size 
by no more than a factor of a few. 

Overall, even if we allow for possible uncertainties dis- 
cussed above, the trend of the ths — relation would not 
be changed. It seems reasonable to suppose that the slope 
change occurs around 1 kpc. 

2.2. Advance speed 

Figure 2 shows the advance speed of hot spots (i'hs)- 
For CSOs, Whs has been directly measured by VLBI ob- 
servations. In the case of CSOs where the hot spot motion 
was measured relative to the counter hot spot (so-called 
the separation rate: ■^sep), we estimated uhs as Usep/2. On 
the other hand, v-hs of FRH sources has not been mea- 
sured directly. We therefore adopted whs of MSOs and 
FRff sources estimated by the synchrotron age constraints 
(Alexander et al.l984; CariUi et al.l991; Liu et al.l992; 
Klein et al. 1995; Mack et al. 1998; Parma et al. 1999; 
Murgia et al. 1998; Schoenmakers et al. 2000; Jamrozy et 
al. 2005), which is estimated by the break frequency of ra- 
dio spectra caused by the synchrotron aging. We adopted 
Whs of MSOs only from the lobe-dominated source (clas- 
sified gtype ah in Murgia et al. 1999). We summarize the 
data of Whs hi Table 2. 




10"! 10" 10^ 

Distance from the core [kpc] 

Fig. 2. — The relation of hot-spot velocity (i^hs) arid hot-spot 
distance from the core (^h)- Pluses (green), stars (blue) and crosses 
(red) represent CSOs, MSOs and FRIIs, respectively. 

As we will show in the foUowings, there are a number 
of biases and uncertainties in the estimation of whs- We 
should conclude that the evolution of whs relative to Ih is 
not clear at this moment. We do not therefore argue the 
favor of the evolutional model based on the observational 
data of Whs • 

Firstly, the apparent absence of CSOs with whs <'~^ 0.1c 
is possibly due to the observational biases. In Figure 3 we 
show the apparent velocity of the hot spot of CSOs with 
respect to redshift of the source. Each parameter is listed 
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in Table 3. Apparent velocity of less than ~ 10 /^as/yr 
has not been detected. This is probably due to a lack of 
angular resolution. Apparent velocity of ^ 10 /J,as/yr cor- 
responds to ~ 0.1c in z > 0.1, so that it is possible to 
conclude that the apparent absence of vns < 0.1c in CSOs 
is due to biases introduced by the detection hmits. 

Secondly, the synchrotron age estimates are more uncer- 
tain compared to the kinematic age estimates. Primary 
uncertainty in the estimation of synchrotron age is the 
magnetic field. Although the minimum-energy condition 
is highly uncertain, the minimum-energy field is usually 
adopted due to the lack of better estimation. Since the 
synchrotron age strongly depends on the magnetic field 



rather than the break frequency (tsyn oc B 



-1.5 -0.5 



un- 



certainty of magnetic field makes it difficult to relate the 
source age. Therefore, whs of MSOs and FRII sources are 
more unreliable than that of CSOs. 
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In this section, we briefly review a dynamical evolution 
model of radio sources (KK06) which traces the dynam- 
ical evolution of advancing hot spots and expanding co- 
coon (e.g., BC89; see also KK05). Hereafter, the mass, 
momentum and kinetic energy of unshocked jet flow are 
conserved. In other words, we do not consider the en- 
trainment effect of the ambient medium. This would be 
justified for the relativistic jet flows (e.g., Scheck et al. 
2002: hereafter S02; Mizuta, Yamada & Takabe 2004). 

The basic equations of the cocoon expansion can be ob- 
tained as follows: 



The equation of motion along the jet axis, i.e., the 
momentum flux of a relativistic jet is balanced to 
the ram pressure of the ambient medium spread 
over the effective cross-section area of the cocoon 
head, A^, 



where Lj, pa aud fns arc the total kinetic energy 
of jets, the mass density of the ambient medium 
and the hot spot velocity, respectively. Here we as- 
sume that Lj is constant in time, and that uhs — ^^h 
where Uh is the advance speed of the cocoon head 
(see KK06 for details). 
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Fig. 3. — Apparent separation velocity between hot spots of 
CSOs. Arrows indicate upper limits (Taylor &; Vermeulen 1997; 
Gugliucci et al. 2005). Concerning the sources in Gugliucci et al. 
(2005), we have not adopted J 1415-1-13 because separation rate is 
measured not between hot spots but jet components. Broken line 
corresponds to the velocity of 0.1c. 

Finally, a direct comparison of -whs between MSOs and 
FRII sources is not acceptable because the estimation 
method of uhs is different between MSOs and FRE sources. 
Basically, whs of FRH sources is calculated by the linear 
fit to the observed synchrotron age distribution across the 
lobe as whs of FRH sources (cf. Fig. 3 in Alexander et al. 
1984). On the other hand, the advance velocities of MSOs 
in Murgia et al. (1999) are estimated in a different way. 
Murgia et al. calculated the break frequency by fitting the 
integrated spectra to the continuous injection model (CI 
model). Then they estimated whs as Ih/tci where tci is 
the synchrotron age estimated from the break frequency of 
the CI model spectrum. The main concern of this method 
is that the integrated spectrum could contain not only the 
lobe emission but also the core and jet emission. To be 
precise, the estimated synchrotron age is the gmeanh par- 
ticle age of these components. The electrons in the core 
and jets, especially in the knots of the jets, are likely to 
be sites of particle acceleration. Thus the estimated syn- 
chrotron age could be underestimated compared to actual 
source age, which leads to the overestimation of whs- 

3. A COEVOLUTION MODEL OF HOT SPOTS AND A 
COCOON 



2. The equation of motion perpendicular to the jet 
axis, that is, the sideways expansion velocity, Vc is 
equal to the shock speed driven by the overpres- 
sured cocoon with the internal pressure, Pc, 

Pc = Pa^c. 

For CSOs, the sideways expansion velocity evalu- 
ated by the outflow velocity of the line-emitting gas 
(e.g., de Veries et al. 1999; Axon et al. 2000; O'Dea 
et al. 2002; Gupta et al. 2005; Holt et al. 2008) 
would be higher than the sound velocity of the am- 
bient medium. Even for FREs, the overpressured 
expansion would be valid (e.g., KK05; Ito et al. 
2008). 



3. The energy conservation in the cocoon, namely all 
of the kinetic energy transported by jets during the 
source age t is is deposited as the cocoon's internal 
pressure, 

PcFe = 2(7e - l)Lji, 

where Vc is the the volume of the cocoon and 
7c — 4/3 is the specific heat ratio of the relativis- 
tic plasma in the cocoon. For simplicity, we do not 
consider the Bremsstrahlung cooling of the cocoon, 
although this would be important in the much ear- 
lier phase of the evolution (Kino et al. 2007). 
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Pa- Also, in order to keep the constant velocity of sideways 
expansion, larger a predicts smaller cocoon pressure Pc be- 
cause of Pc = Pavl- Next we consider the X dependence. 
Larger X leads to larger Ac and Vc, and thus the larger 
Pc is required because of Pc — Pav"^- Thus, the slower vhs 
is needed in order to satisfy the energy conservation in 
the cocoon, that is, PcAcUhs = const. In summary, larger 
(smaller) a and smaller (larger) X leads to acceleration 
(deceleration) of hot spot velocity. 

The aspect ratio of the cocoon, TZ = Ic/hn is also an im- 
portant quantity relevant to the evolution of radio galax- 
ies. The Zh-dependence of the aspect ratio of cocoon is 
then given by 

X(2.5-0.5aj-3 



Fig. 4. — A schematic picture of the evolution of radio sources. 
The aspect ratio of the cocoon is defined as TZ = Ic/lh- 

According to KK06, the physical quantities of hot spots 
and the cocoon can be determined by the mass density of 
the ambient medium, and the cross-section area of the 
cocoon body, Ac (see Figure 4). The mass density profile 
is expressed as 

Paid) oc 

where d is the radial distance from the center of the galaxy. 
The cross-sectional area of the cocoon body is given by 

Ac{t) = nl^ cx t^, 

where Ic = . Vc{t')dt' is the radius of the cocoon body. 
Here tmin is the time from when the two-dimensional (2D) 
phase (^h growth phase) starts. On a free parameter X, 
we can constrain the value of X, in comparison with nu- 
merical simulations. We will discuss this issue later. As- 
suming that ^h/j^HS constant in time, the quantities of 
hot spots (rns , vns and Phs ) can be described in terms of 
the length from the center of the galaxy to the hot spot 
(Zh) as follows: 



(4) 

Finally, in order to check the reliability of their analyti- 
cal model, we compare the results of KK06 compared with 
2D relativistic hydrodynamic simulations in a uniform am- 
bient medium with a — (S02) and in a declining ambi- 
ent medium with a — 1 (Perucho & Marti 2007; hereafter 
PM07) . We summarize the comparisons with 2D hydrody- 
namic simulations in Table 4. These results show that the 
analytic model (KK06) can well agree with the results of 
2D simulations. Thus, it is possible to constrain the free 
parameter X by direct comparison with numerical sim- 
ulations. Following the results of S02 and PM07, we can 
constrain 1.2 < X < 1.4 for the wide range of a (= — 2.5) 
because X does not depend on a significantly (see Table 
4) . As seen in Table 4 we find that the power law index of 
ths is slightly different from the results of KK06 and those 
of numerical simulations (SG2 and PMG7)[3. However, such 
small deviations do not affect the following results. 

4. WHAT CAN WE LEARN FROM OBSERVED Ths — Ih 
RELATION ? 

In this section, we examine the deceleration and accel- 
eration of hot spot velocity by comparing the observed 
ths — Ih relation with the KK06 model. 

4.1. Constant velocity model or constant aspect ratio 

model ? 

Here we test two well-known evolution scenarios of ra- 
dio sources, that is, (i) the constant velocity model (urs = 
const) proposed by Readhead et al. (1996a) and (ii) the 
self-similar model with TZ = const (e.g., Falle 1991; Begel- 
man 1996; Kaiser & Alexander 1997). 

As seen in Fig. 2, one may think this figure seems to sup- 
port the Whs = const model. However, we should keep in 
mind the large uncertainty of uhs (§2.2). Furthermore, the 
deceleration may take place via a strong interaction with 
denser ambient gas in host galaxies, since compact radio 
sources such as CSOs and MSOs would affect a significant 
interaction with the surrounding ambient medium as they 
propagate through it (e.g., Gelderman & Whittle 1994; de 
Vries et al. 1997; Axon et al. 2000; Gupta et al. 2005; Holt 
et al. 2008). At present, it is still under debate whether 

^ The difference from PM07 would come from boundary conditions they adopted. Since they use an open boundary condition in the symmetry 
plane at the jet basis, this allows the leakage of gas from the boundary (see PM07 for more details). Thus, the hot spot pressure drops faster, 
and consequently the growth rate of the hot spot radius is higher. These trends for Phs and rpjs coincide with the deviations between KK06 
and PM07. Concerning the distinction from S02, the effect of relativistic backflow may be a plausible reason considering that the deviation 
appears only for the evolution of the hot spot radius. If the part of jet power is used to fuel the relativistic backflow, its kinetic energy is not 
available to expand the cocoon sideways. As a result, the growth of head and hot spot radius would be faster in S02 than KK06. But the 
deviation (factor 3) is not significant even if we consider the long-term evolution of radio source, e.g., from ~ 10~'' kpc to 10'^ kpc. This would 
imply the relativistic backflow does not have any effect on the cocoon dynamics. 



X(-2 + 0.5a)(c«-2) + 3o 



f HS OC ly^ 



2X(-2 + 0.5a) + 6 



WHS oc 



,X(- 
X(2- 



CX Ih 



0.5a)(a-2) + 4-3Q 



Ph 



X(-2+0.5a) + 3 



CX Zh"" 



(1) 

(2) 
(3) 



and Sr, represent the values of the exponents 



HS OC 

where 5'r, S'v and Sp 
of the hot spot radius, advance velocity and pressure, re- 
spectively. We apply the conditions oiY = X{—2+0.5a) + 
3 > 0, which make the contribution at t^i^ in the integra- 
tion of Ih = . VHs{t')dt' cx — small enough. 
In other words, we do not consider the case that the life- 
times of compact radio sources are much longer than those 
of extended radio sources because this is unrealistic. 

Here, let us consider the a and X dependence on the 
advance velocity of hot spots for fixed physical quantities 
at t = imin- Concerning a dependence in fixed X, larger a 
leads to a weaker deceleration effect on whs due to smaller 
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iifjg = const is really correct. In order to test the validity 
of the Whs — const, model, we predict the required mass 
density profile to explain the observed ths — diagram. 
From eq. (2), the relation of S'v oc 2 - X{2 - 0.5a) = is 
needed for the constant hot spot velocity. By eliminating 
X, a-dependence of ths shows 

''HS OC l^^^. (5) 




0.5 1 1.5 2 2.5 

Slope of ambient mass density (a) 



Pa 



Our model 

(Filled circles) 



a<0.5 




pa 



«=1.5 



VHs=const. model 

(Filled triangles) 

.a -2.6 



a~0.8 



1 kpc 



1 kpc li, 



Fig. 5. — Growth rate of cocoon body {Ac oc t'^) against 
the value of exponents of the ambient mass density profile 
{pa, oc l^")- The blue-shaded region shows the allowed region 
given by the best fit value of the observed rns — relation 
for the CSO-MSO phase {h, < 1 kpc), i.e., = 1.34 ± 0.24. 
The red line corresponds to the allowed region obtained for 
the MSO-FRn phase (Zh > 1 kpc), including the fitting errors, 
i.e., — 0.44 ± 0.08. The black solid line shows iiHS=const 
(Sv ~ 0). The two black dashed lines are relations for the evo- 
lution of hot spot velocity vhs oc l^^ . The light gray shaded 
region represents the range of X-parameter derived from com- 
parisons between KK06 and numerical simulations. The filled 
circles show the solutions for our model. We can constrain on 
he power-law index a{lh < 1 kpc) as < a < 0.5. The filled 
triangles are solutions for the vhs ~ const.. The predicted 
ambient mass density profiles are shown. 

Figure 5 displays the allowed region of the growth rate 
of cocoon body {Ac oc t^) plotted versus a (slope of the 
ambient mass density profile, pa oc ^i7")- The blue-shaded 
region shows the allowed region given by the best fit value 
of the observed ths — relation for the CSO-MSO phase 
(^h < 1 kpc), i.e., S, = 1.34 ± 0.24. The red line corre- 
sponds to the allowed region obtained for the MSO-FRI 
I phase (Zh > 1 kpc), including the fitting errors, i.e.. 
Si = 0.44 ± 0.08. The black solid line shows WHS=const 
(S'v = 0). The two black dashed lines are relations for the 
deceleration and acceleration of fnSj he., S'v = — 1 and 0.3 
(see eq.(2)). The light gray shaded region represents the 
allowed range of X-parameter derived from comparisons 



between the KK06 model and relativistic numerical sim- 
ulations (S02 and PM07). The dark gray shaded region 
shows the forbidden region (see §3). As seen in Fig. 5 
(filled triangles and also the schematic picture), the con- 
stant velocity model can reproduce the observed rns — Ih 
relation only when the slope of the inner part of the density 
profile (< 1 kpc) is steeper than that of the outer part (> 1 
kpc), i.e., a ~ 2.6 (Zh < Ikpc) and a ^ 0.8 (/h > Ikpc). 
However, such a density profile of ambient matter is unre- 
alistic because of the slope of the density profile in many 
clusters of galaxies a « 1.5 ~ 2 (e.g., Trinchieri et al.l986; 
Mulochaey & Zabludoff 1998, Mathews & Brighenti 2003). 
Thus, the vns = const model can be declined. 

Similar to the vus = const model, it is an open ques- 
tion whether the cocoon shape evolves in a self-similar 
way, i.e., TZ = const (De Young 1997; Komissarov & Falle 
1998; O'Dea 1998; Carvalho & O'Dea 2002a, b; PM03). 
Then, we examine whether the TZ = const model can re- 
produce the observed rns — relation. As for the self- 
similar model, the relation of S'r cx X{2.5 — 0.5a) — 3 = 
is needed from eq. (4). By eliminating X, ths can be 
obtained as 

ms oc 4"+^)/^ (6) 

Figure 6 is the same as Fig. 5, but the allowed regions 
for the aspect ratio TZ are plotted. The black sohd line 
denotes TZ = const (S'r, = 0) and the other two black solid 
lines are the solutions of Sr = —0.5 and 0.3. As seen in 
Fig. 6, the self-similar model {TZ = const) cannot explain 
the observed ths — relation because the black solid line 
does not overlap both blue and red regions. Therefore, the 
case of TZ = const can also be ruled out. 




1 1.5 2 2.5 

Slope of ambient mass density (a) 

Fig. 6. — Same as Figure 5, but the black solid line shows 
7^=const (5r = 0) and the two black dashed lines are relations of 
Sr = 0.3 and 0.5 where 7^ oc 

On the whole, two well-known models are not able to 
explain the observed ths — relation, even if we consider 
the possibility that the actual power-law index is steeper 
than Sr > 1.34 ±0.24 for > Ikpc (see §2.1). 

4.2. A new finding: Deceleration and acceleration of hot 
spots velocity 

Here we show an evolutionary track of radio sources 
to be consistent with the observed evolution of ths (see 
Fig. 1). In groups of galaxies and clusters of galaxies 
{Ih > Ikpc), it is well established that the slope of the 
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density profile is a « 1.5 — 2 (e.g., Trinchieri at al.l986; 
Mulochaey & ZabludofF 1998, Mathews & Brighenti 2003)! 
We assume a = 1.5 for /h > Ikpc. Comparing the KK06 
model with the observed rns — relation (Fig. 1), we 
can determine (i) the evolution of the advance speed of 
hot spots, (ii) the slope of mass density distribution for 
Ih < Ikpc and (ill) the evolution of cocoon morphology 
(see Fig. 5 and Fig. 6). 

On the evolution of hot spot velocity, we find that the 
advance speed of the spots and lobes show the decelera- 
tion phase (CSO-MSO phase) and the acceleration phase 
(MSO-FRU phase) as follows (see filled circles showing 
allowed solutions in Fig. 5): In the CSO-MSO phase 
{III < Ikpc), vns decelerates as 

vns oc l^^. (7) 

In the MSO-FRU phase {Ih > Ikpc), vhs slightly acceler- 
ates as 

vm cc (8) 

With respect to the CSO-MSO phase {Ih < Ikpc), the 
flatter density profile (0 < a < 0.5) is predicted in or- 
der to satisfy both the observed ths — hi diagram and 
the constraints from numerical simulations (filled circles 
are allowed solutions in Fig. 5). The predicted exter- 
nal mass density profile is quite similar to a King-profile 
distribution, i.e., pa oc [1 + (r/rc)^]°'^", as indicated by 
X-ray observations of elliptical galaxies (Trinchieri et al. 
1986). Typical core radii in giant elliptical galaxies are 
Tc = 0.5 — Ikpc (Trinchieri et al. 1986). This coincides 
with the observational evidence since most CSO/MSO 
hosts seem to be elliptical galaxies (e.g., O'Dea 1998) al- 
though they show the features of distortions and inter- 
actions. If we consider the mass density distribution we 
predicted, the jets propagate through the denser ambi- 
ent medium during the CSO-MSO phase, compared with 
the MSO-FRU phase. Then, the interaction between the 
jets and the ambient medium is stronger in the CSO-MSO 
phase. This leads to the larger velocity of sideways ex- 
pansion in order to maintain the energy conservation in 
the cocoon. Thus, Ah (or rus) grows faster in the early 
phase of evolution. From the equation of motion along 
the jet axis (sec §3), it is found that the advance speed 
of hot spots {vus) is determined by the linear density of 
the effective working surface, Pa^h, i-e., °^ (Pa^h)"^- 
Hence, the deceleration and the acceleration of fns de- 
pends on the change of the power-law index of ambient 
density profile a in the MSO phase (^ kpc). When pa^h 
increases with as the CSO-MSO phase (pa^h oc l^), the 
hot spot velocity decelerates. On the other hand, the ad- 
vance speed accelerates if pa^h decreases with as the 
MSO-FRn phase (pa^h ^h"'^')- In the present work, we 
impose the idea that the evolution of hot spot velocity is 
mainly determined by the growth of the effective working 
surface, A^. Namely, we did not consider the entrainment 
of the ambient medium in the jets for simplicity (see §3). 
We should stress that it is a new aspect that Ah deter- 
mines the acceleration and deceleration of advance speed 
of hot spots. Note that the study on the entrainment pro- 
cess on jets is also important and complementary for the 
deceleration of relativistic jets within ^ Ikpc (e.g., Bick- 
neU 1984; Laing 1993: De Young 1993; Komissarov 1994; 
Laing 1996; De Young 1997). 



Concerning the evolution of cocoon morphology, we 
show possible solutions (filled circles in Fig. 6). As a 
result, the evolution of a cocoon shape is predicted as fol- 
lows: In the CSO-MSO phase, the shape of a cocoon is 
nearly spherical as 

n oc (9) 

The cocoon shape is elongated as Ih increases in the MSO- 
FRn phase as 

Tl oc 1^°-^. (10) 

Hence, the aspect ratio of cocoon in the MSO phase 
can be close to unity, compared with that in CSOs and 
FRH sources (see Fig. 7). In the CSO-MSO stage, the 
hot spot decelerates with Ih, i.e., t^ns oc l^^, due to 
the strong interaction with the ambient medium. In or- 
der to maintain the energy conservation in the cocoon, 
Pc^c^'HS = const, and the momentum conservation along 
the jet axis, PusAhC = const, the sideways expansion ve- 
locity decreases more slowly compared with the advance 
speed of hot spots, e.g., cx Z^'''^^ for whs cc l^^. Thus, 
the aspect ratio approaches unity since the aspect ratio of 
the cocoon can be described as 7^ « Pc/Phs oc v'c/v^^. In 
the same manner, it is possible to explain the evolution of 
the cocoon shape in the MSO-FRU stage, Vc oc " for 
Vhs oc l^-^. 



CSOs: R - 0.25 




Fig. 7. — A schematic picture of the evolution of a cocoon mor- 
phology derived from our predictions (see Fig. 6). 

Recently, we have noticed a very interesting radio source 
(3C 353) whose radio morphology can be well explained by 
the present model. 3C 353 is an FRU radio galaxy residing 
on the edge of the cluster and has asymmetric radio lobes, 
namely the east lobe is spherical while the west lobe is 
elongated. Thanks to the new XMM-Newton observation, 
it has been found that the east lobe is located in the denser 
region compared with the west lobe (Goodger et al. 2008). 
From the above discussions, the spherical structure of the 
east lobe can be understood due to the stronger interac- 
tion with the ambient medium than the west lobe. If this 
is the case, it is predicted that the advance speed of the 
east lobe and east hot spot is slower than that of the west 
lobe. In future, the advance speed of the east lobe and 
hot spot will approach the sound velocity of the ambient 
medium. Then, the east lobe's type of radio morphology 
will change from FRII into FRI type. Such radio galax- 
ies would be observed as hybrid morphology radio sources 
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(HYMORS), which have clear radio lobes with two dif- 
ferent morphologies (e.g., Gopal-Krishna & Wiita 2000; 
Gawronski et al. 2006). In order to understand the ori- 
gin of HYMORS, i.e., the origin of FRI/FRH dichotomy, 
it would be important to measure the mass density profile 
around HYMORS. 

5. PREDICTION OF FATE OF CSOS: DEAD OR ALIVE ? 

In §4, we found that the observed rjjs ~ relation indi- 
cates that the advance speed decelerates in the CSO-MSO 
stage and accelerates in the MSO-FRII stage. However, 
some CSOs would evolve into some other type of radio 
source like low-power extended radio galaxies (FRI ra- 
dio galaxies), since < 10 kpc scale dying MSOs, which do 
not evolve into FRIIs, have been discovered (Parma et al. 
2007). Thus, we need to further investigate in this section 
the fate of CSOs, or which kinds of CSOs can evolve into 
FRII sources. To this aim, we compare the predicted evolu- 
tion of Whs with the sound velocity of the ambient medium, 
Cs, because the cocoon can expand only when vns > Cg. 
As for the slope of ambient matter density, we assume 
a {Ih < Ikpc) = and a {Ih > lkpc)=1.5 (see §4.2). Cor- 
respondingly, the behavior of whs can be determined as 
'^HS ^ for Zh < 1 kpc and whs oc l^-^ for > 1 kpc (see 
eq. (7) and (8)). The hot ambient-gas temperature, Tg is 
measured in the range of Tg = 5 x 10^ K — 2 x 10'' K, i.e., 
c, = (5fcTg/3mp)i/2 « 7 X 10-^0-1.4 x 10"^^, where is 
the sound velocity of the ambient medium (e.g., Trinchieri 
et al. 1986), k is the Boltzman constant and nip is the 
proton mass. 

Figure 8 shows the evolution of hot spot velocity 
for three initial advance speeds of the ambient medium 
with 'yHs(^h,2D) — 0.01c, 0.1c and 0.5c, where Zh,2D = 
i^Hs(^h,2D)imin IS the distaucc from the core at which the 
2-D phase {A^ growth phase) starts. Here we suppose 
/h,2D = 5pc from Fig. 1. 




10"^ 10"^ 10"' 10° 10' 10^ 10^ 

Distance from the core [kpc] 



Fig. 8. — Our predictions of evolution of i'hSi 'TIS oc l^^ 
for < Ikpc and t>HS oc l^'^ for > Ikpc. The black solid lines 
denote the evolution of hhs to'' ''Hs('h,2D) = 0.01c, 0.1c, and 0.5c 
where ih 2D = 5 pc The black-shaded region represents the evo- 
lutionary path from CSOs into FRUs. The blue-shaded region 
shows the range of sound velocity of the ambient medium, i.e., 
7 X 10-''c < Cs < 1.4 x lO-^c (5 X 10® K < Tg < 2 X lO'^ K) where 
Tg is the temperature of the ambient medium. 

As seen in Fig. 8, we find that CSOs can evolve into FRII 
sources, passing through MSOs for any when WHs(^h,2D) 



is larger than about 0.1c because of WHs(^h) > Cg. On the 
other hand, when UHs(^h,2D) is less than about 0.1c, uhs 
is comparable to the sound velocity of hot spots during 
the CSO-MSO phase {Ih < Ikpc). When uhs equals to Cg, 
the matter and energy of the cocoon mix with the ambient 
medium. After that, the cocoon expands subsonically like 
a sonic boom (Churazov et al. 2000, 2001; Reynolds et 
al. 2001; Brighenti & Mathews 2002; Zanni et al. 2003). 
These results indicate that the advance speed of hot spots 
during the 1 D phase (constant Ah) controls the fate of 
CSOs. Thus, the fate of CSOs can be divided into two 
cases as follows: 

1. The fate of CSO with low WHs(^h,2D) 0.1 c: 
The hot spots velocity is comparable to the sound 
velocity of ambient matter at the MSO stage (< 
1 kpc) . Consequently, the ambient gas may be 
heated up since the relativistic plasma within the 
cocoon leaks into the ambient medium. The cocoon 
morphology is distorted in the MSO phase. After 
that, radio sources would expand with the sound 
velocity. Thus, the evolutionary sequence appears 
as follows: 

CSO MSO FRI. 

2. The fate of CSO with high WHs(^h,2D) 0.1 c: 
The CSOs can evolve into FRII sources, passing 
through MSOs whose cocoon shape is spherical. 
Namely, the evolutionary sequence is as follows: 

CSO ^ MSO ^ FRH. 

The above predicted cocoon morphology of MSOs may 
coincide with observations in which the fraction of dis- 
torted morphology of MSOs is higher than that of CSOs 
(Saikia et al. 1995; O'Dea 1998; Dallacasa et al. 2002a, 
b; Parma et al. 2007). Moreover, with respect to case (1), 
Drake et al. (2004) observed radio-weak CSOs from the 
IRAS survey, then they found that radio-weak CSOs are 
young (< 10^ yr) and show slightly irregular morphology 
radio sources (also see Kunert-Bajraszewska et al. 2005). 
Therefore, such radio-weak CSOs may be young counter- 
parts of ~ 1 kpc scale low-power compact radio sources 
(Giroletti et al. 2007). Such low-power CSOs might be 
observed as the dying MSOs (e.g., Parma et al. 2007). As 
another possibility, low-power CSOs may be the progeni- 
tors of FRIs because the radio lobes expand with w Cg, and 
then the projected linear size of radio sources can reach 
~ 30kpc(= 10"^c X 10* yr) if the typical age of FRIs is 
w 10* yr (Parma et al 1999). If this is the case, the origin 
of the FRI/FRH dichotomy could be related to the dif- 
ference in initial advance speed of hot spots uhs,id which 
depend on Lj and pa (see §3). 

Moreover, these findings impact on AGN feedback asso- 
ciated with AGN jets. The influence of AGN activity on 
the evolution of the host galaxy is now a hot and timely 
topic. It has been suggested that the AGN bubbles quench 
star formation and regulate the SMBH growth, and conse- 
quently the local BH-to-bulge relation can be well repro- 
duced (e.g.. Silk & Rees 1998; King 2003). However, we 
should recall that their conclusions rely on the dynamical 
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evolution of AGN bubbles with the constant shock veloc- 
ity {vhs = const), assuming a singular isothermal hot- 
ambient distribution (pa oc r~^). However, as we men- 
tioned the constant velocity model does not match with 
our predictions derived from the actual evolution of young 
radio sources. Therefore, we must reconsider the condi- 
tions that AGN feedback due to AGN jets works effec- 
tively. This is left for our future work. 

Finally it is of value to mention what we should do as 
future observations, to reveal the dynamical evolution of 
radio sources as follows: 

1. We need to evaluate the accurate v^s of MSOs in 
future VLBI observations, to test the deceleration 
and acceleration of hot spot velocity. At the same 
time, it is important to estimate the advance ve- 
locity of hot spots for FRUs by kinematic studies 
because most vas for FRlIs have been derived by 
synchrotron-aging work only. 

2. In order to confirm the evolution of the cocoon 
shape we predict here (see Fig. 7), it is worth sys- 
tematically measuring the aspect ratio of cocoons 
for CSOs, MSOs and FRUs. But the emission of 
the cocoon body is very faint in general due to the 
synchrotron cooling, and thus it would be hard to 
assess the clear shape of a cocoon (i.e., Ac) by the 
present facility. The sensitivity and low frequency 
capability of SKA can be used for direct measure- 
ment of cocoon morphology, since the cooling effect 
is smaller at lower frequency. Instead of the direct 
measurement of the cocoon shape, it is also impor- 
tant to measure the morphology of radio lobes, i.e., 
Alf^/lh (see Fig. 4) as far as Ac oc Ah is a reason- 
able assumption. 

3. In the present paper, we impose the idea that the 
initial advance velocity of hot spots, t'Hs(^h,2D), 
may be interpreted by the differences in total kinetic 
energy. However, we cannot rule out the possibility 
that the ambient mass density controls i'Hs(^h,2D)- 
Thus, it is crucial to systematically explore whether 
there is a significant difference in ambient mass den- 
sity for CSOs, MSOs and FRIIs. For this purpose, 
the search for cold gas via molecular gas and 21 cm 
HI is very useful using future facilities with high 
spatial resolution such as ALMA and SKA. At the 
same time, the study for diffuse hot gas will be also 
important using future X-ray satellites with high 
sensitivity and high spatial resolution such as XUES 
and Constellation-X. 

4. In order to reveal the end of the ID phase of hot 
spot evolutions, it is necessary to discover and elu- 
cidate the candidates for ultra compact symmetric 
radio sources (e.g., Ih < 5pc) such as High Fre- 
quency Peakers (HFFs) that have been intensively 
observed and a growing number of the sources are 
being detected (e.g., Dallacasa et al. 2000; Tinti et 
al. 2005; Orienti et al. 2007), using high resolu- 
tion VLBI such as VSOP-2 (see also Snellen 2008). 
This is very important not only to predict the fate of 



CSOs as we mentioned, but also to reveal the phys- 
ical conditions for the transitions from ID phase to 
2D phase. 

5. It is worth exploring whether radio- weak CSOs, in- 
cluding radio-quiet AGNs (Seyferts and radio-quiet 
quasars) follow the same rns — diagram. The 
present model predicts that the observed ths — Ih 
relation of bright radio sources holds radio-weak 
sources as far as the AGN jets are composed of the 
relativistic plasma. So far, it has been very difficult 
to measure the hot spot properties (e.g., size and 
velocity) because of the faintness of low-powered 
jet sources. However, future low-frequency radio 
telescopes (e.g., LOFAR and SKA) are ideal tools 
for understanding the dynamical evolution of radio- 
weak sources from the point of view of sensitivity 
and the low-frequency capability. 

6. SUMMARY 

We have investigated the physical relationship between 

CSOs, MSOs, and FRIIs, together with the coevolution 
model of hot spots and a cocoon (KK06) and the obser- 
vation reflecting the dynamical evolution of radio sources. 
For a larger sample of radio sources than previous work 
(JOG; PM03), we find that the evolution of the hot spot 
radius clearly changes at /h ^ Ikpc, even after we carefully 
exclude the observational bias. Based on these observa- 
tional characteristics of hot spots, we investigate the evo- 
lutionary track of extragalactic radio galaxies. The main 
results of the present work are summarized as follows: 

(i) We test two well-known evolution models of radio 
sources, i.e., "the constant velocity" model and "the con- 
stant aspect ratio" model. Since the mass density pro- 
file of the ambient matter shows unrealistic behavior, i.e., 
the derived power law index {a « 0.8) is shallower than 
a = 1.5 — 2 in groups of galaxies and clusters of galaxies, 
the constant velocity model can be ruled out. In the case 
of the self-similar model, it turns out that the constant 
aspect ratio model cannot reproduce the observed ths ^ 
relation. Therefore, the case of constant aspect ratio can 
also be declined. 

(ii) We investigate the evolution of hot spot velocity to 
reproduce the observed rns — relation. As a result, we 
find that the advance speed of hot spots and lobes strongly 
decelerate when the jets pass through the ambient medium 
in host galaxies (i.e., the CSO-MSO phase), while the jets 
slightly accelerate outside host galaxies (i.e, the MSO-FRI 
I phase). Recent observation shows that a possible corre- 
lation between the outflow velocity of the ionized gas and 
Ih might indicate a possible deceleration of the jets as they 
pass through the host galaxy (Labiano 2008). This trend 
seems to be consistent with the preferred model, although 
the available data is not enough to compare with our pre- 
diction. Furthermore, we can constrain the allowed value 
of a as < a < 0.5 for Zh < 1 kpc. The main origin of the 
deceleration and acceleration of hot spots is the change 
in power-law index of ambient density distribution at ~ 
1 kpc. It is worth emphasizing again that the advance 
speed of hot spots is determined by the growth of the co- 
coon head, which is a new aspect of the deceleration and 
acceleration of relativistic jets. 
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(iii) It is also found that the cocoon shape of MSOs be- 
comes spherical or disrupted, while the cocoon morphology 
of CSOs and FRIIs is elongated. The radio image of MSOs 
seems to be consistent with our prediction (Dallacasa et 
al. 2002a, b). However, the emission of the cocoon is gen- 
erally very faint due to the synchrotron cooling and thus it 
is quite difficult to observe it using present facilities. Thus, 
it is of value to systematically measure the aspect ratio (or 
the morphology of radio lobes) for CSOs, MSOs, and FRIIs 
using future low- frequency radio telescopes (e.g., LOFAR 
and SKA). 

(iv) We predict the fate of CSOs, comparing the advance 
speed of hot spots with the sound velocity of the ambi- 
ent medium. Only CSOs whose initial advance speed is 
higher than about 0.1c can be progenitors of powerful FRI 
Is. It would be in agreement with the observation which 
the fraction of FRIIs is less than that of CSOs. Moreover, 
the fate of CSOs would be closely related to the FRI/FRI 



I dichotomy and AGN feedback which we will discuss in 
our up-coming paper. 
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Kawakatu, Nagai and Kino 



Table 1 

Physical parameters of hot spots 



see http://yso.mtk. nao.ac.jp/'"kawakatu/Kawakatu08.pcif 



Table 2 
Advance speed of Hot spots 



see http://yso.mtk. nao.ac.jp/'^kawakatu/KawakatuOS.pdf 



Table 3 

Apparent separation rate (wapp) between hot spots 



Source 


z 


Vapp [MasJ 


Reference 


0108+388 


0.699 


9.27+1.21 


(1) 


0710+439 


0.518 


17+1 


(2) 


1031+567 


0.4597 


14.6+4.8 


(3) 


1245+676 


0.1071 


34.9+1.8 


(4) 


1943+546 


0.263 


25+2 


(2) 


2352+495 


0.238 


32.7+11 


(3) 


OQ208 


0.0766 


32+20 


(5) 


CTD93 


0.473 


23+7.7 


(6) 


3C84^ 


0.0183 


270 + 70 


(7) 


Jllll+1955 


0.299 


< 10 


(8) 


J1414+4554 


0.19 


< 14 


(8) 


J1734+0926 


0.61 


< 8 


(8) 


1946+708 


0.101 


< 60 


(9) 


J1934-638^ 


0.183 


< 30 


(10) 



^"Uapp is measured between the core and hot spot. 

References. — (1) Owsianik, Conway & Polatidis 1998, (2) Polatidis & Conway 2003, (3) Taylor et al. 2000, (4) Marecki et al. 2003, (5) 
Wang et al. 2003, (6) Nagai et al. 2006, (7) Asada et al. 2006, (8) GugUucci et al. 2005, (9) Taylor & Vermeulen 1997, (10) Tzioumis et al. 
1989 



Table 4 

Comparison with 2D relativistic hydrodynamic simulations 





WHS 


ms 




-Phs 






(i) A uniform ambient medium 


a = 0)" 




S02 


7-0.55 
'h 


;0.45 
'h 




KK06 


7-0.56 
'h 


;0.55 
'h 




'h 


PM07 
KK06 


7-0.11 
7-0.11 


(ii) A declining ambient medium 

;0.78 
'h 

;0.67 
'h 


{a = 1)"^ 


pl.4 



NOTE.- "The case (i) corresponds to KK06 model with X — 1.2, while the case (ii) corresponds to KK06 model with 
X = 1.4 



